Analog gravity from electrodynamics in non-linear media 
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Abstract 

Working with electrodynamics in the geometrical optics approximation we derive the expression 
representing an effectively curved geometry which guides the propagation of electromagnetic waves 
in material media whose physical properties depend on an external electric field. The issue of 
birefringence is addressed, and the trajectory of the extraordinary ray is explicitly worked out. 
Quite general curves are obtained for the path of the light ray by suitably setting an electric field. 

PACS numbers: 04.20.-q, ll.lO.-z 
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I. INTRODUCTION 



The theory of general relativity lead us to the con- 
ception of a generally curved geometry of the spacetime 
whose source encompass all matter and energy content 
of the universe. The comparative study between the 
kinematic aspects of general relativity and other kinds 
of interactions has recently been called as analog models 
for gravitation W. fndeed, it have extensively been ex- 
amined, not only with respect to electrodynamics (non- 
linear |§-0 as well as in moving dielectric fluids [|-12|) 
but also for acoustic perturbations |13|-[l5|. Long ago, 
Gordon pq ] shown that the refractive index of a medium 
can be reinterpreted by means of an effective metric 
Scordor. — v'^'^ ~ (1 ~ eii)V^V'^ , where V^ represents the 
velocity 4-vector of an arbitrary observer. The paths of 
light rays in non-linear electrodynamics can be obtained 
in terms of the geodesic equations of an effective geom- 
etry representing the non-linearities 0,0] ■ These works 
assumed electrodynamics as governed by a quite general 
non-linear Lagrangian which is a function of the two in- 
variants of the Maxwell field. 

Inside material media. Maxwell equations must be sup- 
plemented by constitutive laws that relate the electro- 
magnetic excitations and the field strengths by means of 
quantities characterizing each medium the waves prop- 
agate into. In this context, the electromagnetic inter- 
action was geometrized for the case of linear constitu- 
tive relations ||l^;|l^. More recently, specific non-linear 
constitutive relations were also geometrized liy,EG]. The 
above results lead us to conclude that electromagnetic 
waves propagate inside material media as if immersed in 
a curved spacetime. This fact allows one to make an 
analogy between wave propagation in non-linear media 
and gravitational phenomena. 

This work explicitly presents, in the context of geo- 
metrical optics, the construction of an effective geometry 
which describes small perturbations of the electromag- 
netic field inside substances whose electric susceptibility 
depend on the electric field. A description of birefrin- 



gence [Q observed for this class of materials is presented, 
and the problem of the trajectory of the light ray is ex- 
plicitly presented in a similar way as the bending of light 
due to gravitational interaction. 

We work in Minkowski spacetime employing a 
Cartesian coordinate system. The background met- 
ric will be represented by 77^^, which is defined by 
diag{+\, —1, —1, —1). Units are such that c = 1. 



II. FIELD EQUATIONS 

i^'"' and P^'^ are the tensors representing the total 
electromagnetic field, which are expressed in terms of the 
strengths and the excitations of the electric and magnetic 
fields as 






(1) 
(2) 



where the Levi-Civita tensor introduced is defined such 
that 77°i23 = -f 1. 

In general the properties of the media are determined 
by the tensors e'^p and /i"/3 which relate the electromag- 
netic excitation and the field strength by the generalized 
constitutive laws, 






(3) 

(4) 



In the absence of sources Maxwell's theory can be sum- 
marized by the equations 



P^'\^ = 
h'" , = 



(5) 
(6) 



together with the constitutive laws (ph and (0), where we 
have introduced the dual electromagnetic field tensor 



p^.. ^ _^^^^^p^P^ 



(7) 



Eq. (Q) could be alternatively obtained by the Bianchi 
identity. Now, rewriting the above field equations in 
terms of the electric and magnetic fields it results 






(8) 
(9) 



Additionally, the electromagnetic excitation is related 
to the field strength by means of the constitutive rela- 
tions (H) and (Q), whose derivatives with respect to the 
coordinates can be presented as 















E^H^\r (10) 

H'^H^'^r- (11) 



III. GEOMETRICAL OPTICS 

In order to determine the propagation of waves asso- 
ciated to the electromagnetic field, we will consider the 



method of field discontinuities [£1| . We define a surface 
of discontinuity E by z(a;^) = 0. Whenever S is a global 
surface, it divides the spacetime in two distinct regions 
[/" for z{x>^) < 0, and t/+ for z{x'^) > 0. The disconti- 
nuity of an arbitrary function /(x") on E is given by 



[/(^")], 



lim \f(P+) 



fiP- 



(12) 



with P+, P^ and P belonging to C/+, U^ and E, respec- 
tively. The electric and magnetic fields are continuous 
when crossing the surface E. However, their derivatives 
behave as 



[E^,,U = e^K,; [H^.A 



hf'K,, 



(13) 



where e^ and h^ represent the discontinuities of the fields 
on the surface S and 



K, 



as 



(14) 



is the wave vector. Applying these conditions to the field 
equations (pi) and (0), we obtain 



J 



ef^pe^ ■ 






9e" 



dE'' 



Ef^Kr, 



dey 



Ef^K^hf" 



^e^' 



dE^ 



^E^e^ 



9£^ 



dEt' 



H^Kr. 






dH'^ 



E^h'' {KV) + v^^c^pK^V^h" 



[f^'ph^ + %^HP^.- + ^^'^") ^KV) - r,^'\,K.V-e^ = 



(15) 
(16) 
(17) 
(18) 



where we have defined {KV) = K^^Vp^. Eqs. (p^-(p^) 
came from the zeroth component of Eqs. (Bll^CTi ^^'^ 
can be obtained from Eqs. (pT[)-(p^) upon contraction 
with K^ whenever {KV) ^ 0. Then, it follows that the 
coupled system ([l7|)-([l8|) completely describe the prop- 
agation of light rays. 



J 



For the cases Sap — Sap{E^,H^^) and fia/S = fJ-{Vap 
V°'V^) with /i a constant, Eqs. (17)-(|lq) reduces to 



/i" 



KKV) 
and the eigen-vector equation 



T^^P-PKpV^ep 



^(^^)'(^''/' + ^^" 



{KV)^7^^^"pEPKrV„ - K^Kp + {KV)V^Kp + KH^^ - {KV)H^; 



(19) 



(20) 



r 



The Fresnel equation represents non-trivial solutions of 
Eq. (pfl) and is given by the determinant of the term in 
brackets. The solution of the Fresnel equation for linear 
constitutive relations is derived in reference |[1§1. 



Now, let us analyze non-linear media whose physical 
properties are influenced by an external electric field as 



.M/3 ^ e(^>)(^M/3 _ v'^V'^) - a{E^)EfEf^. 



(21) 



Each particular medium is characterized by the param- problem reduces to 
eters a and e. Considering Eq. (pi|), the eigen- vector 

J 

I [K^ - {KVf + n{KVf{e + aE^)] S^^ - K^^Kp + {KV)V''Vp - £>" (2aEp 



de 



where we denoted E 



'^-dE0 

-E°'Ea- Expanding the polarization vector in a suitable basis 
e^ = aE^ + hH^ + cK^^ + dV^ 



E'^)^e» = ,22) 



(23) 



one gets 



K' 



Wvf-T^'^ "'^'' + ^"(11^ + ^'^)]+' h^'"' + "■' {^ 



-2aE^KR + K'^ I - 



de 



E' 



.-] 



= 



\dEi3 dEP ) 

b [K^ - {KVf + n{KVf{€ + aE"^)] = 
aE^^K^, + bH^'K^, + c{KVf [l - ^{e + uE"^)] + d{KV) = 
a{KV)E^'K^ + b{KV)mK^ + c{KV)K'^ + d [K^ + iJL{KVf{€ + aE^)] = 0. 



E' 



da \ 
dE0) 



(24) 

(25) 
(26) 
(27) 



There are two distinct solutions of the above system. 
First, by fixing b ^ 0, it follows that a = c = d = 0, 
and we obtain the light cone condition 



K"^ = {KVf 



n 

' E 



(28) 



where il = E{e + aE^). Conversely, assuming 6 = 
and a ^ 0, it follows that c = [E Ef^Kf3/nn{KVf] a and 
d = —[EE'^Ki3/fifl{KV)] a, and results in the light cone 
condition 



K^ = {KVf[l 



E^ dn\ EE'^K^K'^ d fVl 



M- 



E dE0 



n 



dEP\E 



(29) 



Therefore, there are different light-cone conditions for 
each polarization eigenvectors e^ . Corresponding to Eq. 
( Pq ) one has the polarization mode 



bH^ 



traveling with phase velocity 



< = 



E 



(30) 



(31) 



Moreover, Eqs. (|24|), (|2§)-(g7D for this mode require the 
following conditions: H^Kfs = 0, a{E^Hx){E'^Kg)^ 0, 
and HP{d/dE'^){n/E) = 0. Similarly, from Eq. (|2|) the 
polarization mode and its phase velocity are given by 



^jlEP 



^Mj<vf_n 

E^Kx E 
^ 1 



on 



El^ + Kl^ - {KV)V'^ 



n ^ dE<^ \E 



(32) 
(33) 



where we introduced the space-like unit propagation vec- 
tor fc" = (ry"'3 - V°'V^)K0/^{KVf - K^ . Contrarily 
to the velocity v^, the velocity v_ depends on the di- 
rection of the wave propagation. Both these velocities 
coincide when the wave propagates along the direction 
of the electric field Ef^ {i.e., P = ±E^ jE). The ap- 
pearance of two different modes with different velocities 
points out to birefringence phenomena ||2^,|2^, in this 
case induced by the electric field. The modes (pOJ ) and 
(P2) correspond to the ordinary and the extraordinary 
rays, respectively. The corresponding refraction indexes 
are given by n± = l/w±. 



IV. ANALOG GRAVITY 



3|) can be restated in the more appealing 



Eqs. (|2 
form 



g^^'^K^K, = (34) 

where we defined the symmetric tensors 



3^ = 77^"^ - ( 1 - /^-^ 1 Vf'V 



g"- = v""- 



E 



ep on \ 

' E dEP 



^ ^ ^^\ /(M^-) 



2VtdEP \E 



(35) 



(36) 



in which a'^'^'^-' = a'^'^ + a'^^ for an arbitrary tensor a^''. 
The inverse symmetric tensor g^^ is defined in such way 
that g'^^ga. = <5^ 

Differentiation of Eq. (p4|) with respect to x'^ yields 



2K,^^K^g'''' + K^K,g''\ 



0. 



(37) 



As the vector K'^ is a gradient (04) it follows that 
-f^/^.A = Kx^^. Moreover, one has the identity 



.g^^A 



-g'^'^gt^" (Sctcr, A + ffoA, <t - 5(tA, 

Thus, Eq. (|37| ) can be presented as 

g^''K^[Kx,.-r''.xK^]^0 
where we defined the quantity 



r^crA = ■7^g^°'{gaa,\ + ga\,a - ga\,a)- 



(38) 
(39) 

(40) 



Eqs. ( |35| ) and (|6|) depend on the total field E, which 
encompasses both the external field and the one associ- 
ated with the wave itself. Thus, the standard treatment 



of birefringence [p2 23 1 by setting a frequency-dependent 
e"^ p{uj) is included in the present formalism as a partic- 
ular case as d/duj = {dE"^ /duj)d/dE'^. We will hereafter 
adopt the point of view that the wave fields are much 
smaller than their external counterparts. Thus, the fields 
appearing in g^'^ are independent of the light ray. 

Therefore one recognizes that, for the propagation vec- 
tor K^, the tensor g'^" is a metric tensor. Indeed, K^, is a 
light-like (or null) vector with respect to this metric ten- 
sor, Eq. (34), and also satisfies the geodesic equation (|39| ) 
in terms of the Christoffcl symbols r"^,y associated with 
this metric, Eq. (|o|). We thus refer to K^, as a geodesic 
null vector with respect to the effective metric tensor g^''. 
The two terms in brackets in Eq. (|3^) correspond to the 
covariant derivative of K^, and will henceforth be denoted 
as K^-\. The effective geometry g^^'^ allows us to write 
K^ — g^^Ki,, for which the geodesic equation can be 
written in the simpler form 



Rf^.^K 



A 







(41) 



which avoids explicit dependence on the metric tensor. 



V. TRAJECTORIES 

The trajectories of light-rays are the geodesic lines 
which solve the geodesic equation (Ba). For the sake 



of simplicity, we will consider V^ = 5'^. Only the cal- 
culations for the extraordinary ray will be presented, 
as the corresponding ones for the ordinary ray follow 
exactly the same reasoning (being much simpler, how- 
ever). In order to fit it best problems with axial sym- 
metry, standard cylindrical coordinates a;^ = {t, r, 9, z) 
will be adopted. For a configuration of electric field 
£:m = (0, E) with E = {Er, Eg, E^) = E^[r, 9)z, and 
assuming dVl/dEr = = dil/dEg, the metric associated 
with the extraordinary ray reduces to 



^fj-i^ 



f^ E dE^ 



E, dn 

E dE^ 






-1 
















1 
-pi 














o 



dn 

n dE^ 



(42) 



Thus, g^^ = diag[fi{r, 9), -1, -r^, -/2(r, 9)], where we 
introduce the notation 



fiir,0)^ M 



E^ dfl 
's'dEl 



h{r, 9) = 



E^ dn 
TldEl 



(43) 
(44) 



We will also denote the components of the propagation 
vector as K'^ = (dt, dr, d9, dz)/dT, where t is an affine 
parameter along the integral curves of the vector K'^ , 
and [t, r, 9, z) represent the coordinates of the points 
belonging to these lines. The geodesic lines for this case 
are given by Eq. (|l]) as 



dH 

dT^ 


1 


9/1 dxP dt 
dxP dr dr 






(45) 


d^r 

dT^ 


1 

^ 2 


dh fdzV dh fdtV 
dr \dT ) dr \dT J 


deV 

\dTj 


(46) 


d^9 
dr"^ 


1 

2^2 


"a/2 fdzV dh (dtV 
39 \dT) d9 \dT) 


2dr de 
r dr dr 


(47) 


d^z 
dr^ 

The 


1 9/2 dxP dz 
f2dxf> dr dr' 

solutions of Eqs. (4|) and (48) are 




(48) 




dt A 




(49) 








dz B 






(50) 



where A and B are arbitrary constants. Inserting these 
relations in Eqs. (E^)-(E^) one gets 



(51) 
(52) 



d^r fd9Y 


1 d r 


A^ 


B^^ 


dr"^ "^ \dTj 


^ 2dr 


h 


/2J 


d'^e 2 dr de 

dr^ r dr dr 


1 3 
2r^de 


p2 
./i 


52 

/2 



For dr/dr ^Q^ dO/dr equations (|lD-(|52|) lead to 



VI. CONCLUSION 



dr\' 
dr) 



dev 

d^J 



h 



/2 



D 



(53) 



where D is an arbitrary constant. Introducing Eqs. ([49|)- 
(pO|) and d53) in the Ught-cone constraint 



ds^ = g^^dx'^dx" = g^^K''K"dT^ 







(54) 



it follows that D = 0. Thus, Eq. (|53D yields dO/dr in 
terms of r which is given by 



dr2 



ld_ 
r dr 



A^ 



B^ 
h 



(55) 



Cylindrical symmetry 



For the simple case /i = 
from Eq. (H) that 

d6 
'db- 



/i (r) and /2 = /2 {r) it follows 



C_ 



with an arbitrary constant C, from which Eq. 



(56) 
33|) yields 

(57) 



The trajectories of the light ray is then found from the 
ratios {dz / dT) / [dr / dr) and {dO / dr) / {dr / dr) , from which 
one obtains z = z{r) and d = 6{r) as 




±B 



e = ±c 



dr 



A^ B^ 
dr 



C^ 



A^ B^ 



c2 



(58) 



(59) 



Eqs. ([49D-(pOD, ( |53| ) and ( ^5D encompass in this way a 
complete description for the path of the light ray with 
arbitrary ^(E^) = E[e{Ez) + a{Ez)E'^] as a function 
of Ez{r, 9) in terms of the arbitrary constants A, B, C. 
Thus, quite general curves can be obtained for the path 
of the light ray by suitably setting the electric field E^ 
at the given medium. From u* = {dx^ / dr) / {dt / dr) it 
follows that 



E 



^^^cTE: 



[2f,A^J 



f dVt 



\dE, 



(60) 



which is to be compared with Eq. (|33|). Birefringence 
then occurs for this model whenever \dQ,/dEx\ / |ri/i?2;|, 
which can be ascribed either to a non zero a parameter 
or to a varying e{E) in equation (|2l|), or even to both of 
them. 



Maxwell linear equations together with the constitu- 
tive relations for non-linear materials lead to an effec- 
tive non-linear theory of electrodynamics. The problem 
of the propagation of electromagnetic waves was here 
dealt with under the approximation of geometrical op- 
tics. The path of light rays thus belong to the char- 
acteristic surfaces which depend on both the external 
fields and the propagation fields. For the case of a con- 
stant scalar magnetic permeability a Fresnel problem is 
obtained |20|. The propagating modes and their veloc- 
ities were obtained for the class of electric permittiv- 
ity e^'^{E). Weak waves were then described as light 
rays propagating along null geodesies of an effectively 
deformed geometry which depends on the external elec- 
tric field, and birefringence phenomena were found to oc- 
cur if de^'^ /dE^ ^ 0. Explicit calculations for the path 
and velocity of the rays were presented for the particular 
e/^" = e{E)[r]t"' - V^'V] - a{E)E''E'' case for cylindri- 
cally symmetric media, illustrating the usefulness of the 
geometrical approach. The results of this investigation 
may be useful to test kinematic aspects of general rela- 
tivity in the laboratory. Indeed there are many works in 
the recent literature dealing with such theme, not only in 
the context of light propagation in dielectric media; see 
for instance the review in Ref. [[131 . 



ACKNOWLEDGMENTS 



This work was partially supported by the Brazil- 
ian agencies Conselho Nacional de Desenvolvimento 
Cientifico e Tecnologico (CNPq) and Fundagdo de Am- 
paro a Pesquisa do Estado de Minas Gerais (FAPEMIG). 



[1] Workshop on Analog models of general relativity, 
www.cbpf.br/~bscg/analog, Rio de Janeiro, Brazil, Oc- 
tober 2000. 

[2] S. A. Gutierrez, A. L. Dudley and J. F. Plebanski, J. 
Math. Phys. 22, 2835 (1981). 

[3] J. S. Heyl and L. Hernquist, Phys. Rev. D59, 045005 
(1999). 

[4] W. Dittrich, and H. Gies, Phys. Rev. D 58, 025004 
(1998); ibid, Phys. Lett. B 431, 420 (1998). 

[5] M. Novello, V. A. De Lorenci, J. M. Salim and R. Klip- 
pert, Phys. Rev. D 61, 45001 (2000). 



[6] 

[7] 
[8] 
[9] 

[10] 

[11] 
[12] 

[13] 

[14] 
[15] 
[16] 
[17] 



V. A. De Lorenci, R. Klippert, M. Novello and J. M. 

Salim, Phys. Lett. B 482, 134 (2000). 

S. Liberati, S. Sonego and M. Visser, Phys. Rev. D63, 

085003 (2001). 

U. Leonhardt and P. Piwnicki, Phys. Rev. A60, 4301 

(1999). 

U. Leonhardt and P. Piwnicki, Phys. Rev. Lett. 84, 822 

(2000); M. Visser Phys. Rev. Lett. 85, 5252 (2000); U. 

Leonhardt and P. Piwnicki, Phys. Rev. Lett. 85, 5253 

(2000). 

C. Barcelo, S. Liberati and M. Visser, Clas. Quantum 

Grav. 18, 1137 (2001). 



L Brevik and G. Haines, gr-qc/0106045 . 

C. Barcelo, S. Liberati and M. Visser, gr-qc/0104001. 

W. G. Unruh, Phys. Rev. Lett. 46, 1351 (1981); Phys. 

Rev. D51, 2827 (1995). 

T. Jacobson Phys. Rev. D48, 728 (1993) 

M. Visser, Clas. Quantum Grav. 15, 1767 (1998). 

W. Gordon, Ann. der Phys., 72, 421 (1923). 

Y. N. Obukhov and F. W. Hehl, Phys. Lett. B 458, 466 



(2000); M. Schonberg, Revista Brasileira de Fisica 1, 91 

(1971). 
[18] Y. N. Obukhov, T. Fukui and G. F. Rubilar, Phys. Rev. 

D 62, 4405 (2000); F. W. Hehl, Y. N. Obukhov and G. 

F. Rubilar, a:r-qc/9911096 A. Gross and G. F. Rubilar, 

gr-qc/0102022 . 
[19] M. Novello and J. M. Salim, Phys. Rev. D63, 083511 

(2001). 
[20] V. A. De Lorenci and M. A. Souza, Phys. Lett. B512, 

417 (2001). 
[21] J. Hadamard, in Legons sur la propagation des ondes et 

les equations de I'hydrodynamique (Ed. Hermann, Paris, 

1903); J. Plebanski, in Lectures on non-linear electro- 
dynamics (Nordita, Copenhagen, 1968); G. Boillat, J. 

Math. Phys. 11, 941 (1970); 
[22] L. Landau and E. Lifchitz, Electrodynamique des milieux 

continus (Ed. Mir, Moscou, 1969). 
[23] M. Born and E. Wolf, Principles of optics (Cambridge, 

6*'' edition, 1980). 



